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We calculate correlation functions of the energy-momentum tensor in the vicin- 
ity of the deconfinement phase transition of (3-|-l)-dimensional SU(2) gauge theory 
and discuss their critical behavior in the vicinity of the second order deconfinement 
transition. We show that correlation functions of the trace of the energy momen- 
tum tensor diverge uniformly at the critical point in proportion to the specific heat 
singularity. Correlation functions of the pressure, on the other hand, stay finite at 
the critical point. We discuss the consequences of these findings for the analysis 
of transport coefficients, in particular the bulk viscosity, in the vicinity of a second 
order phase transition point. 

I. INTRODUCTION 

Experimental evidence for rapid thermalization of the dense matter created in heavy ion 
collisions at RHIC [1] provided strong motivation for a new look at the non-perturbative 
structure of the quark gluon plasma (QGP). The rapid thermalization of the QGP has been 
related to the smallness of the shear viscosity in the plasma phase [2] , which may come close 
to the theoretical lower bound established for the ratio of shear viscosity to entropy density 
[3]. This gave rise to the interpretation of the quark gluon plasma above but close to the 
transition temperature as a strongly interacting medium that has properties of an almost 
perfect liquid. 

These experimental findings also renewed the interest in determining transport properties 
of gauge theories through the calculation of correlation functions of the energy-momentum 
tensor on the lattice. This is by far not a straightforward calculation and requires a careful 
analysis of the low frequency structure of the spectral representation of these correlation 
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functions [4]. Nonetheless progress has been made through high statistics calculations in 
SU(3) gauge theories [5-7]. 

Recently it has been argued that close to the transition from low temperature hadronic 
matter to the plasma phase of QCD bulk viscosity might play a much more important role 
than shear viscosity [8]; bulk viscosity has been related to a temperature derivative of the 
trace of the energy-momentum tensor. The latter does diverge at a second order critical point 
[9] and thus would force bulk viscosity to diverge too. This argument is based on a sum rule 
derived for correlation functions of the energy- momentum tensor in Minkowski space [8]. 
However, as noted recently [10], in Euclidean space the correlation function of the trace of 
the energy-momentum tensor picks up an additional constant contribution corresponding to 
an exact 5-function at zero frequency in the spectral representation of this correlator. This 
contribution, in fact, cancels the leading singular behavior of the correlation function and 
thus makes the relation between transport coefficients and properties of bulk thermodynamic 
observables more subtle. 

In the following we will provide evidence for the presence of a contribution to the corre- 
lation function of the trace of energy-momentum tensor that is constant in Euclidean time 
but temperature dependent. We show that the temperature dependence of this contribution 
scales like the specific heat and thus dominates the singular behavior of this correlation 
function in the vicinity of T^.. These findings also explain the inconsistency observed in the 
parametric dependence of bulk viscosity derived in high temperature perturbation theory 
and through the sum rule analysis [11]. 

The singular behavior of bulk viscosity in the vicinity of a critical point has long been 
known in statistical physics [12, 13]. In particular, the divergence of bulk viscosity at the 
critical point of the liquid gas transition has been analyzed extensively [14]. In this case the 
divergence of bulk viscosity [Q has been related to the critical behavior of the 3-dimensional 
Ising model [13, 15, 16] , which is in the universality class for the critical point of the liquid-gas 
phase transition. It has been argued that the divergence of C is strong and, in fact, almost 
quadratic in the inverse reduced temperature, t = \T — T^jTc- The singular behavior, 
~ t"^^"*"", with a, V denoting static critical exponents of the 3-d Ising model and z being 
a dynamical exponent characterizing the equilibration of density fluctuations, goes along 
with a strong divergence of the relaxation time for density fluctuations, t/j. Their ratio, 
C/tr ~ however, is proportional to the inverse of the specific heat and thus vanishes 
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slowly at the critical point. 

While numerous numerical studies have confirmed that the critical behavior of bulk ther- 
modynamics in (ri + l)-dimensional gauge theories is indeed controlled by c?-dimensional 
universality classes corresponding to the global symmetry of the relevant order parameter, 
little is known about the relation between dynamic properties and the related critical expo- 
nents [17]. The SU(2) gauge theory with its second order deconfinement phase transition 
seems to be an ideal model to explore critical behavior of dynamical properties, e.g. trans- 
port coefficients. The SU(2) gauge theory at finite temperature incorporates all the basic 
features one expects to be relevant for deconfinement in QCD. Moreover, it has a second 
order phase transition which belongs to the universality class of the 3-dimensional Ising 
model. It thus will allow to analyze in how far dynamic universal properties known for the 
Ising universality class are relevant for transport properties in a quantum field theory. It 
also may give insight into transport properties in the vicinity of the chiral critical point in 
QCD that may exist at non-zero baryon number density and also would belong to the Ising 
universality class [18]. In the vicinity of the SU(2) deconfinement phase transition a lot of 
information exists about the critical behavior of bulk thermodynamics that can be related 
to the critical behavior of correlation functions of the energy- momentum tensor. This will 
help to gain experience with lattice calculations of the bulk viscosity in pure gauge theories. 

In this paper we will establish relations between the critical behavior of bulk thermody- 
namic observables in the vicinity of the SU(2) deconfinement transition and properties of 
correlation functions of the energy-momentum tensor. In particular, we will analyze the Eu- 
clidean correlation functions of the trace of the energy-momentum tensor and compare it to 
critical behavior of the specific heat. We present evidence that the strength of the divergence 
of this correlation function at is independent of Euclidean time and is controlled by the 
universal structure of the divergence of the specific heat at a critical point belonging to the 
3-dimensional Ising universality class. In the next section we review the calculation of bulk 
thermodynamic observables in SU(2) gauge theory [19] and present some new results on the 
diagonal components of the energy-momentum tensor, i.e. energy density and pressure, in 
the vicinity of the deconfinement transition temperature. In Section 111 we introduce local 
operators for energy density and pressure and discuss basic properties of their finite temper- 
ature Euclidean time correlation functions at vanishing momentum. Section IV is devoted 
to a discussion of critical behavior extracted from correlation functions of the trace anomaly 
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and presents a calculation of the pressure-pressure correlator. We analyze consequences for 
the determination of bulk viscosity from these correlation functions in Section V. Finally, 
we conclude in Section VI. 

II. THERMODYNAMICS OF SU(2) GAUGE THEORIES ON THE LATTICE 

As preparation for the analysis of correlation functions of diagonal elements of the energy- 
momentum tensor we want to discuss here the calculation of the relevant bulk thermody- 
namic observables, energy density (e) and pressure (P), on the lattice. In order to analyze 
Euclidean space-time correlation functions of e, P or the trace of the energy-momentum 
tensor, = e — 3P, one obviously needs to use an approach to bulk thermodynamics 
that allows to define local operators for energy density and pressure. We will call this ap- 
proach, which in fact was the basis for the original formulation of bulk thermodynamics on 
the lattice [21], the differential formahsm. It is in contrast to the so-called integral method 
which nowadays is most commonly used in lattice calculations to extract energy density and 
pressure from the trace anomaly [22] . The differential formahsm has not been used recently 
in studies of bulk thermodynamic observables as it is more involved and requires additional 
non-perturbative calculations of derivatives of the bare gauge couplings with respect to the 
anisotropy parameter, ^ = Qa/ar, that controls the spatial (a = a^) and temporal (ttr) 
lattice spacings [23] . Having control over the anisotropy is needed to derive thermodynamic 
quantities as function of temperature, T — and volume, V — {Nf^Qa-)^ through ap- 

propriate partial derivatives with respect to a,- and [21], although numerical calculations 
will generally be performed on isotropic lattice — 1). In particular, for the study of the 
singular behavior of the energy density in the vicinity of the SU(2) deconfinement transition 
this approach has been used successfully [19]. In this context also the necessary derivatives 
of the gauge coupling with respect to ^ had been analyzed in detail. We will follow this 
analysis here closely. 
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A. The differential formalism 

We start from the partition function for a SU{Nc) gauge theory at finite temperature, 
which we write in standard lattice notation [24], 

Z{T,V)^ /^dt^^.Me-"^^■^^ (1) 

x,n 

with 

O AT 

Here x — {x, X4) labels the discrete set of space-time points on a four-dimensional hypercubic 
lattice of size N^N^.. 

In the standard Wilson discretization scheme the local gauge action in the /iz/-hyperplane 

is expressed in terms of plaquette variables P^i,{x) = -^'R.eTTUx,^Ux+ef,,uUl^g^ ^Ul with 
link variables Ux,fx G SU{Nc). Wc write the action in terms of contributions arising only 
from entirely space-like plaquettes {P^) and space-time like plaquettes {Pt), 

X i>j=l,2,3 

=l^Y.Pr{x) , Prix) = IY1 • 
X i=l,2,3 

We also introduce Pq — [Pa- + Pt)/'^ as the average of space- like and time-like plaquettes 
evaluated on a symmetric (T ~ 0) lattice, i.e. for N^- — N^. 

Prom the logarithm of the partition function one obtains pressure and energy density as 
derivatives with respect to and a^, respectively [21], 

^ = N,N^ {{2g-' - (C, - Cr)) (Pa - Pr) - 3 (C, + Cr) (2Po - (Pa + Pr))) , (4) 

^ = 3N,N^ {{2g-^ - (c, - c,)) (P, - P,) + (c, + c,) (2Po - (P, + P,))) . (5) 

Here and Cr are given in terms of derivatives of the gauge coupling, g~^{a, {), with respect 
to the anisotropy at = 1 [21, 23]. Their sum is constraint by the /3-function, 



Big-') = a 



dg 



-2 



= -2 (C, + Cr) . (6) 

^=1 



da 

With these relations for e and P it is straightforward to obtain also the trace of the 
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energy-momentum tensor and the entropy density, 

6B{g-^)N,N^ [P, + Pr- 2Po] , (7) 



rp4 rp4 

s _ e + p 



= 47V^iV,^ ( [2g-' - (c, - c.)] {P^ - P,)) . (8) 



We note that these particular combinations of energy density and pressure are proportional 
to simple sums and differences of space-like and time-like plaquette expectation values. 
They only require common multiplicative renormalizations for P^r and P^. Extracting e/T^ 
or P/T^ directly, however, involves different renormalization factors for space-like and time- 
like plaquette expectation values. For this reason it is straightforward to analyze correlation 
functions of local densities of the entropy or trace anomaly. We note that in the vicinity of the 
SU(2) deconfinement transition the dominant singular behavior of ©^'^/T^ as well as s/T^ 
is controlled by the energy density. We will show later that in the vicinity of Tc the pressure 
is an order of magnitude smaller than the energy density. When one combines 
and s/T^ to extract the pressure it thus requires good knowledge of the renormalization 
factors Co- and in order to eliminate the dominant energy contributions and arrive at a 
proper definition for a local operator for the pressure. This is a prerequisite for the analysis 
of pressure-pressure correlation functions and, as will become clear later, this becomes of 
importance for studies of transport coefficients, e.g. the bulk viscosity, in the vicinity of Tc. 



B. Critical energy density and pressure 

As an application of the differential formalism reviewed above we have performed a finite 
size scaling analysis of the critical energy density and the pressure at the deconfinement 
transition temperature of the SU(2) gauge theory. This extends earlier studies of the critical 
energy density performed in Ref. [19] to larger lattices and yields, for the first time, a 
determination of the pressure at the critical point. 

The deconfinement transition of the SU(2) gauge theory in (3+l)-dimensions is in the 
same universality class as the 3-dimensional Ising model. Its critical behavior is well un- 
derstood and has been confirmed in various finite size scaling studies. On a lattice of size 
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FIG. 1: Finite size scaling analysis of bulk thermodynamic observables on lattices of size A',- = 4, 6 
and 8. The Curves show fits that have been performed for N-r/Na- < 1/6. 



N^Nt energy density and pressure at the critical point are expected to scale as 



\ T4 



/ Nr,N^ 



C 



+ ap 



{l-a)/u 



(9) 
(10) 



Here a = 0.110(1) and u = 0.6301(4) denote the critical exponents of the 3- dimensional 
Ising universality class, controlling the divergence of the specific heat and the correlation 
length, respectively^. 

We have analyzed the scaling behavior of energy density and pressure on lattices with 
temporal extent Nr = 4, 6 and 8 for various spatial lattice sizes; covering the range iVo- G 
[16, 96] for A^^ = 4, and A^"^ G [24, 96] for = 6 and A^^ G [32, 96] for A^^ = 8. For these 
values of Nr the infinite volume critical couplings are known quite precisely [25-27]. We have 
listed them in Table I. At these values of the gauge coupling also the anisotropy coefficients 
Co- and Cr have been determined previously [19]. 

We have calculated pressure, energy density, entropy density and trace anomaly on lat- 
tices of temporal extent A^.^ = 4, 6 and 8 for various spatial lattice sizes at the critical 
couplings given in Table I. Data for the vacuum subtraction needed for all observables but 



Note that the leading finite volume correction that arises from the singular contribution to the pressure is 
inverse proportional to the volume as it is the case also for the leading corrections to regular terms. This 
is due to a hyper-scaling relation, 2 — a = diy with d = 3. 
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the entropy density have been obtained through simulations on a 32^ lattice. The num- 
ber of configurations analyzed ranges from 1 • 10^ on our smaller lattices to 4 • 10^ on the 
large lattices. For our calculations on = 4 and 6 lattices we have used a standard heat- 
bath/ovcr-rclaxation algorithm. For calculations on the larger A^,- = 8 lattices it became 
possible to use also a two-level algorithm [20] which turned out to be more efficient for the 
analysis of correlation functions. Autocorrelation times in the transition region range from 
0{1) on the small lattices to about 40 on the largest lattices. 

In Fig. 1 we show results for the volume dependence of various bulk thermodynamic 
observables. It clearly can be seen that the scaling behaviour of the pressure is different from 
that of the other three observables. Moreover, it is apparent that the pressure is more than 
an order of magnitude smaller than any of the observables that contain contributions of the 
energy density. This shows that the finite size scaling behavior of {e+P)/T^ and (e — 3P)/r^ 
is dominated by the non-analytic behavior of the energy density and thus is expected to be 
controlled by the same ansatz as for the energy density given in Eq. 10. We also note that 
cut-off effects seem to be small already for lattices with temporal extent Nt>6 and that the 
volume dependence starts being controlled by the scaling variable V^^^T = N^j/Nr- This 
may be expected as bulk thermodynamics in the vicinity of is dominated by low frequency 
modes for which the underlying lattice cut-off becomes unimportant. 

To check consistency with the scaling behavior of the Ising universahty class we have 
fitted the data for the energy on the Nr — ^ lattice for all N^/Nr > 4 to the ansatz given in 
Eq. 10 with a free scahng exponent ^ — {1 — a)/v. This fit yields A = 1.41(6) which is in 
good agreement with the value in the 3-dimensional Ising universality class, A = 1.412(1) 
[29]. Fits to the finite volume corrections for the pressure yield an exponent d = 3.0 ± 1.4 
in good agreement with the expected hyper-scaling relation d = 3. 

To extract the asymptotic, infinite volume values for the critical energy density and the 
pressure we then used the ansatze given in Eqs. 9 and 10 with fixed scaling exponents 1.412 
and 3, respectively. Results for e{%)/T^ and P{Tc)/T^ for N^/N^ = 12 as well as infinite 
volume extrapolated values are summarized in Table I; results for the critical energy density 
are consistent with those of Ref. [21] but have a statistical error that is an order of magnitude 
smaller. Comparison of the Nt- = 8 results with those for N^- = 6 suggest that the remaining 
cut-off dependence is small. 



9 





Pc 


/ f-r—t \ t r-r—tA 




/ p-w—i \ t r-r—iA 










UNr 




CO 


4 


2.29895(10) 


0.3303(9) 


0.0240(2) 


0.2872(5) 


0.0242(3) 


6 


2.4265(30) 


0.3015(27) 


0.0136(8) 


0.2722(31) 


0.0135(8) 


8 


2.5115(40) 


0.3158(42) 


0.0107(26) 




0.0107(15) 



TABLE I: Critical energy density and pressure calculated on lattices with temporal A^,- on lattices 
with aspect ratio N^/N^ = 12 and infinite volume extrapolated results. The second column gives 
the values of the gauge couplings at which these calculations have been performed. The error 
estimates are based on fits with the ansatze given in Eqs. 10 and 9 and does not include errors 
that arise from uncertainties in the couplings c^, Cr as well as the critical coupling fic- 

III. CORRELATION FUNCTIONS OF THE ENERGY-MOMENTUM TENSOR 

Within the differential framework for bulk thermodynamics on the lattice, reviewed in 
the previous section, it is straightforward to define local operators for the energy density 
and pressure, as well as the entropy density and the trace of the energy momentum tensor. 
Using Pt{x) and Pa{x) as introduced in Eq. 3 one obtains from Eqs. 4, 5 and Eqs. 7, 8 local 
operators for pressure, energy and entropy density as well as the trace anomaly, e.g. 

e^'^(f , X4) = 6N,B{g-^) [P,{x, X4) + Prix, X4) - 2Po] (11) 

for the trace anomaly and similar for other observables. The local plaquette variables, 
Pa{x,X4), Pt{x,X4) have been defined in Eq. 3. We will denote this particular choice of 
local expressions for plaquette variables in the following as discretization scheme 1. However, 
when considering correlation functions in Euclidean time, one may want to account for the 
effect that time-like plaquettes connect two different time-slices in a lattice and should be 
connected to a Euclidean time that is displaced from that of space-like plaquettes. This 
ambiguity arises because gauge fields are introduced on links of the lattice and the field 
strength tensor thus is spread at least over a domain of the size of an elementary plaquette. 
This effects introduce additional O(a^) discretization errors, i.e. errors that are of the same 
order as the discretization errors that exist anyhow in the definition of the Euclidean action 
or the energy-momentum tensor used in our calculations. They, of course, will disappear in 
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the continuum limit. Nonetheless, comparing different discretization schemes will allow to 
estimate the magnitude of systematic C(a^) effects. We therefore have introduced two other 
discretization schemes which symmetrize the contribution of space-like plaquettes relative 
to a time-like plaquette or vice verse, 

scheme 2: Pij(x) \ {Pij{x) Pij{x 64)) , 
scheme 3: PQj{x) — > \ {P4j(x) + Poj{x — 64)) . 

Note that there is no need for a local definition for Pq which enters the definition of G'*'* {x) . 
Pq is evaluated at zero temperature and is a constant, identical at all space-time points. We 
furthermore introduce zero-momentum projected operators and their fluctuations, 

" X 

Y{t) = Y{t)-{—Y,Y{t)). (12) 

Correlation functions are then obtained as thermal averages over products of fluctuation 
operators, i.e. we consider connected correlation functions, 

^''yl^^^ = Nl{X{T)Ym . (13) 

Due to the periodic boundary conditions in Euclidean time points in the hyperplane at 
T = are identical to those at r = A^,-- In the following we also will discuss properties of 
the midpoint-subtracted correlation functions, 

"^^^y^^^ = Nl {{X{t)Y{0)) - (X(l/2T)y(0))) . (14) 

We used the three discretization schemes introduced above for diagonal components of 
the energy-momentum tensor to analyze the cut-off dependence of correlation functions. In 
Fig. 2(left) we show results for correlation functions of the trace of the energy momentum 
tensor. Gee and the energy, G^e- These correlators have been evaluated on lattices with 
the largest temporal extent, A^^ = 8, used in this study, and at the largest Euclidean 
time separation possible at finite temperature, i.e. the midpoint of the temporal direction 
r = 1/2T. It is apparent from these figures that cut-off effects in Gee ^-re only of the order 
of 10%; within the current errors different discretization schemes yield consistent results on 
the Nr — S lattices. For Gee, however, discretization errors are large even on the largest 
lattices. 
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FIG. 2: Cut-off dependence in different discretization schemes of the energy momentum tensor. 
The left hand part of the figure shows results for the trace-trace and energy-energy correlation 
functions for different values of the lattice cut-off aT = l/N^ evaluated on lattices with spatial 
extent N^- = 12Nt-. The right hand part of the figure shows results for the energy-energy and 
energy-pressure correlation functions versus Euclidean time, tT. 

In the right hand part of Fig. 2 we show the energy-energy and energy-pressure correlation 
functions as function of Euclidean time, rT, evaluated at the critical temperature. As has 
been pointed out in Ref. [10] these correlation functions will, in the continuum limit, be 
independent of Euclidean time as the energy operator is time independent, 

GMr,T) = T'-^{Xm. (15) 

At finite lattice spacing, however, this is not yet the case. It is apparent from Fig. 2(right) 
that the influence of a finite lattice cut-off is more severe at shorter distances, and correlation 
functions G^x thus show a significant dependence on Euclidean time. As a consequence, 
the simple relation between midpoint-subtracted correlation functions of B^^ and pressure- 
pressure correlators, 

AGee{r,T) = AG,,(r,T) + 9AGpp{t,T) - 6AG,p(r,T) (16) 
= 9AGpp{T,T) + 0{a^) , (17) 

which holds in the continuum limit, is violated by O(a^) corrections at finite lattice spacing. 
The left hand side of Eq. 16 as well as the various components of the right hand side are 
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FIG. 3: Midpoint-subtracted correlation functions on lattices with temporal extent = 8 versus 
temperature. Shown are various contributions to the correlation function of the trace anomaly 
(Eq. 16). 

shown in Fig. 3 for tT = 1/4 and A^,- = 8. In the continuum limit the differences AG^t:{T, T) 
and AG^p{r,T) should vanish and AG'0e(T, T) should be equal to 9AGpp(r, T). This is 
clearly not the case for finite values of the lattice cut-off. 



IV. CORRELATION FUNCTIONS AND CRITICAL BEHAVIOR 

Fortunately, the cut-off dependence discussed in the previous section is not crucial for the 
analysis of critical behavior as it arises from short- distance effects. Close to Tc, however, the 
correlation length is large and thermal effects are not very sensitive to the underlying cut- 
off. The cut-off dependence thus is part of the smooth regular background that contributes 
to Gxy{'^,T). This is apparent from the weak temperature dependence of the differences 
'^Gxy{'^,T) shown in Fig. 3. Moreover, we find that correlators calculated in the three 
discretization schemes show a similar temperature dependence and differ, at fixed tT only 
by an almost temperature independent shift. 



A. Correlation function of Qt^t^ 



We will discuss here the temperature dependence of correlation functions of Q^^{x, r) at 
vanishing momentum, i.e. we calculate Gqq{t,T) /T^ as introduced in Eq. 13. This corre- 
lation function is the sum of three contributions corresponding to energy-energy, pressure- 
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pressure and energy-pressure correlations, respectively, 

Geeir, T) = G,,{t, T) + 9G'pp(t, T) - QG,p{t, T) . (18) 

As discussed above in the continuum limit the entire r-dependence of Gee is expected to 
arise from the pressure-pressure correlations, while close to Tc the dominant temperature 
dependence will arise from the energy-energy correlator, which through the fluctuation- 
dissipation theorem, is proportional to the specific heat, 

Geeir, T) ^ de ^ Cy 

Close to the deconfinement transition we therefore expect that Gee(''";^) will show critical 
behavior that coincides with that of the specific heat in a 3-dimensional Ising model, 

T-T, 



Geeir, T) cy 

J^5 J-S =■= 



T-T, 



l + 5± 



+ ... for T^T^. (20) 



Here a is the specific heat critical exponent and uj is the correction to scaling exponent 
which characterizes the leading non-analytic scaling correction. Like a also the exponent 
UJ is universal. For the SU{2) gauge theory in (3-|-l)-dimensions, which we are analyzing 
here, the relevant universality class is that of the 3-dimensional Ising model. In this case 
a — 0.110(1) and u; — 0.53(3) [29] which also has been verified in lattice calculations [19]. 

At non-zero lattice spacing the direct relation between correlation functions involving the 
energy operator and temperature derivatives of any observable is violated by cut-off effects. 
Nonetheless we expect that these are small in the vicinity of a second order phase transition. 
One thus may expect that at least the singular behavior of correlation functions that involve 
correlations with the energy operator will be independent of Euclidean time. Will analyze 
the critical behavior of Geeir, T) in the next two subsections and show explicitly that the 
leading singular behavior indeed is independent of Euclidean time. 



B. Finite size scaling of Gee at the critical point 

We have calculated the correlation function Geeir, T) close to the deconfinement tran- 
sition point of the SU(2) gauge theory. In most of our simulation we use lattices of size 
N^Nr with Nr — 4. This gives us information on the correlation function at two non-zero 
values of Euclidean time, i.e. at tT = 1/4 and at the midpoint tT — 1/2. Most of our 
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a 


V 


uj/v 


9A 


0.110(1) 


0.6301(4) 


0.84(4) 


0.54(1) 



TABLE II: Specific heat exponent (a), correlation length exponent (z^), correction-to-scaling expo- 
nent u and the universal ratio of amplitudes qa = ^+/^- of the 3-dimensional Ising model (from 
[29]). 

simulations have been performed at temperatures close to the phase transition where the 
correlation length becomes large. This required calculations on large spatial lattices in order 
to eliminate finite volume effects. We used spatial lattice sizes with aspect ratios N^/Nt 
varying from 8 (32^ x 4 lattices) up to values as large as 32 (128^ x 4 lattices). This large 
aspect ratio made it possible to perform calculations at temperatures in the proximity of T^. 
without being affected by large finite volume effects. We checked that for all temperature 
values, except for calculations performed directly at T^, the thermodynamic limit has been 
reached within our numerical accuracy. To reach the statistical accuracy needed for our 
scaling tests close to the critical point, where fiuctuations are large, we have generated a 
large number of gauge field configurations ranging from 1 • 10^ on our smaller lattices to 
4 • 10^ on the large lattices. 

As mentioned in Section TLB the location of the critical point is quite well known for 
lattices with temporal extent A^,- = 4 [25]. The relation between gauge coupling and temper- 
ature has been determined in the vicinity of j3c non-perturbatively from calculations of the 
transition temperature on lattices of temporal size 4 < A^ < 16 [19] . We use these results 
to determine the reduced temperature t. 

In various studies of the critical behavior of bulk thermodynamics of the (3-1-1)- 
dimensional SU(2) gauge theory it has been established that this gauge theory belongs to the 
universality class of the 3-dimensional Ising model. In our analysis of the critical behavior 
of correlation functions of the energy-momentum tensor we use established results on the 
critical exponents a, v and uj as well as the universal ratio of the specific heat amplitudes, 
Aj^/A^. These parameters are summarized in Table II. 

In Fig. 4(left) we show results for the correlation function Gqq{t, T)/T^ calculated at the 
midpoint, tT — 1/2, on lattices with temporal extent N^- — A and for various spatial lattice 
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FIG. 4: The correlation function of the trace of the energy-momentum tensor, Gqq{t,T)/T^ at 
distance tT = 1/2 calculated on lattices with temporal extent A^^^ = 4 and several spatial lattice 
sizes No- (left). The curves show spline interpolations. The right hand part of the figure shows 
the volume dependence of Gqq{t,T)/T^ at the critical point T = Tc- The curves show finite size 
scaling fits as explained in the text. The hardly visible open squares on top of the tT =1/2 data 
are the data obtained for tT = 1/4 shifted by ACo- = 18.0675. 

sizes. As can be seen at all temperatures except, of course, at Tc the dependence on the 
spatial volume is well under control. The peak at the pseudo-critical point close to Tc rises 
only slowly with increasing size of the volume as the specific heat exponent a is quite small. 
The right hand part of Fig. 4 shows the volume dependence of the correlation function 
evaluated at Tc on lattices with temporal extent Nr = 4 at the two non-zero distances 
accessible in this case. At Tc we have fitted the volume dependence of G@0{t,Tc)/T^ to a 
finite size scaling ansatz, which includes the leading singular behavior of the specific heat 
and a correction-to-scaling term, 

Gee(r, Tc)/Tc^ = A^N^/'^ (l + B^N^'^/'') + . (21) 

Here Ag., B^, C„ are fit parameters, which all might depend on Euclidean time. We per- 
formed fits for the two non-zero Euclidean time separations reachable on lattices with tem- 
poral extent A'^ = 4, i.e. tT = 1/4 and 1/2. In the fits we make use of the known Ising 
exponents a, v and lo which are given in Table II. This ansatz gives excellent finite size 
scaling fits for both data sets shown in Fig. 4(right). This confirms that the critical be- 
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tT 


A+ 


B+ 




C 


D 








free 


-1 1 A 

1/4 


16.3(2.1) 


-0.59(19) 


-1.9(1.4) 


21.5(5.7) 


-150(53) 


9.2(1.2) 


-2.4(1.0) 


26.3(2.7) 


fit 


1/2 


16.4(2.2) 


-0.76(18) 


-2.2(1.5) 


3.7(3.4) 


-95(32) 


9.1(1.2) 


-2.4(0.9) 


8.4(2.9) 


combined 
fit 


1/4 
1/2 


16.5(1.3) 


-0.74(10) 


-2.12(86) 


21.8(2.1) 
3.8(1.9) 


-143(20) 
-87(18) 


9.15(73) 


-2.39(59) 


26.4(1.9) 
8.3(1.7) 



TABLE III: Fit parameters for fits to Gqq{t,T)/T^ using Eq. 22 and for the volume dependence 
of Gqq{t,Tc)/T^. The last row gives parameters of combined fits to the data for tT = 1/4 and 
1/2. 

havior is controlled by the exponents of the 3-dimensional Ising model. When performing 
fits of the data sets at distance rT = 1/4 and 1/2 separately, we observe that the ampli- 
tudes Arj and B^r agree within errors. This also is apparent from the fact that the data at 
tT — 1/4 can be shifted on top of the data at tT = 1/2 by a independent constant 
(open squares in Fig. 4(right)). We then performed simultaneous fits of both data sets, still 
allowing for different amplitudes at both distances. Prom this combined fit we find that 
A^(rT = 1/2)/^^(tT = 1/4) = 1.01±0.19 and B^(rT = 1/2)/B^{tT = 1/4) = 1.01±0.58. 
The fit results^ are summarized in Table III, which also includes the coefficients extracted 
from an equally good combined fit of both data sets obtained requiring the same coefficients 
A+, B^ and B^ for the two Euclidean distances. 

The analysis presented above estabhshes that the correlation function of the trace of the 
energy-momentum tensor shows the expected universal singular structure of the specific heat 
in a 3-dimensional Ising model. In the vicinity of Tc the singular contributions to Gqq{t, T) 
are found to be independent of Euclidean time. We have performed corresponding analyzes 
for Nt- = Q and 8 which are consistent with these findings, however statistical errors are 
larger in these cases. 

The analysis presented above confirms that the correlation function Gee contains a con- 
stant contribution that arises from low frequency modes that lead to the singular behavior 



^ Note that the coefficients of the correction-to-scahng term, B±, are negative. Although these amplitudes 
are non-universal it has been noted previously that for the 3-dimensional Ising model on various lattice 
geometries this correction term comes out to be negative [30]. 
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FIG. 5: The correlation function of the trace of the energy-momentum tensor, Gqq{t,T)/T^ at 
distances tT =1/4 and 1/2 (left). The curve shows the combined fit of both data sets performed 
with the ansatz given in Eq. 22. The dashed straight lines give the regular part of the fit. The 
right hand part of this figure shows the difference AGee = G@q{1/4:T,T) — Gqq{1/2T,T). 

of the energy-energy correlation function at Tc. In the continuum limit this constant will 
reflect the r-independence of the correlators G^e and G^p. 



C. Critical behavior of G&&{t,T) in the vicinity of Tc 

In the previous subsection we established through a flnite size scaling analysis at Tc that 
the singular terms in GeQ{T,T) are independent of Euclidean time. We will show here that 
this is also the case at temperatures in the vicinity of at which the universal scaling 
behavior allows to separate the long-distance singular behavior of Ge@{T,T) from regular 
contributions. 

In Fig. 5 (left) we show inflnite volume extrapolated results for Ge@{T,T)/T^ at two dis- 
tances, tT = 1/4 and 1/2. Here we left out the data at Tc. At both distances Gee{T,T)/T^ 
increases rapidly when T approaches Tc, above Tc the correlation function changes by a 
factor 2 in the temperature interval [1.01 Tc, 1.1 TJ. 

In the vicinity of Tc we expect that the data are well described by the scaling ansatz, 

Gee(r, T)/T'' = ^±^-"(1 + B±r) + G + Dt, (22) 



where A^, B±, G, D are free parameters, which again all may depend on Euclidean time; 
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— QaA^ and a, u and qa are the known Ising exponents given in Table II. We again have 
performed fits for the two data sets at rT = 1/2 and 1/4 separately as well as simultaneously. 
The ansatz given in Eq. 22 gives good fits in the interval T/Tc G [0.94, 1.05]. As expected 
the parameters C extracted for tT = 1/2 and 1/4, are consistent with the constant terms 
obtained from the finite size scaling analysis at Tg. 

Again we find that within errors the fit parameters and B± are independent of r; a 
combined fit of GeQ{l/AT,T) /T^ and Gee{l/2T,T)/T^ with common amplitudes A+ and 
B± also gives a x^/dof = 1.1. The resulting fit parameters are also given in Table III. 

We thus conclude that in the vicinity of the critical point of the SU(2) gauge theory the 
temperature dependence of the correlation function of the energy- momentum tensor reflects 
the singular behavior of the specific heat of 3-dimensional Ising model. More striking may be 
the fact that the amplitudes for the singular terms are to a good approximation independent 
of T. However, as discussed above this is consistent with the expected r-independence of 
the energy-energy correlation function [10]. 

In order to eliminate the leading singular behavior from Gqq (t, T) it thus suffices to 
consider differences of the correlation function evaluated for different Euclidean time sepa- 
rations, e.g. 

AGee ^ G'ee(l/4r,r) _ Gee{l/2T,T) 

Fig. 5(right) shows that this difference indeed stays finite at T^. However, as shown in Fig. 2, 
at non-zero value of the lattice spacing AGee does not yet coincide with 9AGpp as it will 
in the continuum limit. 

We finally show in Fig. 6 the temperature dependence of the different correlation func- 
tions, Gxx{t.iT), contributing to Gee- As can be seen the entire strong temperature de- 
pendence of Gee in the vicinity of Tc indeed arises from G^e- All other correlators show 
only a weak temperature dependence in the vicinity of Tg. In particular, we find that the 
pressure-pressure correlation function stays finite at Tc and varies little in its vicinity. As 
will be discussed in the next section this has immediate consequences for the determination 
of the bulk viscosity from correlation functions of the energy-momentum tensor. 
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FIG. 6: Correlation functions Gxx{^/'^T^T) contributing to the correlator of the trace anomaly. 
Shown are results obtained from calculations on lattices with A',- = 4 in the infinite volume limit. 

V. BULK VISCOSITY AND RELAXATION TIME SCALE 

Correlation functions of the energy-momentum tensor are widely used to gain informa- 
tion on transport coefficients. The spectral representation of these Euclidean-time correla- 
tion functions and their relation to real-time Kubo formulas allow, in principle, to extract 
transport coefficients from the zero frequency limit of spectral functions [4]. We want to 
discuss here some consequences for the calculation of transport coefficients, in particular the 
bulk viscosity, that can be drawn from the analysis of correlation functions presented in the 
previous sections. 

The bulk viscosity can be extracted from the spectral representation of pressure-pressure 
correlation functions [4], 

poo 

Gpp{t,T)= / du;ppp{iJ,T)K{r,u;,T), (24) 
Jo 

where 

cosh (uj(t — ^)) , ^ 

^ ' ' ' smh(^) ^ ' 

and ppp{uj,T) is the spectral density at vanishing momentum. Making use of the definition 
of transport coefficients via the Kubo-formulas it is readily seen that the bulk viscosity ( is 
related to the low frequency limit of this spectral function [28], 



(26) 
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Prom Eq. 18 it is obvious that this can also be obtained from a spectral analysis of the 
correlator of the trace of the energy-momentum tensor. However, as we have verified explic- 
itly in the previous sections this correlation function receives a constant contribution from 

correlation functions containing the energy operator. This constant term is reflected in zero 
frequency contribution to pee, ^-c a 5-function at vanishing frequencies, 

/■oo 

Gee{T,T)= / duj pee{uJ,T)K{T,u,T) , (27) 

with 

peeico, T) = 9ppp(a;, T) + Tdrie - ^P)u5{u) . (28) 

As the bulk viscosity is obtained from the spectral function in the limit of vanishing frequency 
one should eliminate the contribution arising from pee(O). This can, of course, easily be 
done by analyzing the subtracted correlation function 1\Gqq{t,T) introduced in Eq. 23. 
Prom Pig. 5 (right) we know that this correlation function will stay finite at Tc. Its spectral 
representation is given by 

/■°° . cosh (w(t - 75^.)) - 1 , , 

Jo smD.[^) 

The low frequency part of the spectral function, ppp{uj,T), is expected to be linear, the 
slope being proportional to the bulk viscosity. The rise of the spectral functions ends at 
a characteristic frequency cuq, which is identified as a characteristic relaxation time scale, 
Tp = I/cuq, for density fluctuations. Of course, cuq will be temperature dependent and, in 
fact, is expected to vanish at a second order phase transition point where fluctuations grow 
and relaxation times will diverge. 

We may split the spectral function ppp into a low frequency part and a large frequency 
part, ppp{u!,T) = f{u!,u!o) + P>{i^,T). The low frequency behavior of the spectral function 
may be modeled by a Breit-Wigner ansatz, 

/(a;,^o) = §C(T)(^-5X7^ ■ (30) 

With this we can write the low frequency contribution to the pressure correlator, Gpp{tT, T), 
as 

.low(.rr.^ 9^,. r ^ox cosh(a;oa;(r - 1/2T)) - 1 



71 Jq I + x'^ smh(a;oa;/2T) 
= 9rC(T)u;o(r) (l - ^ In [2 - 2 cos(27rrT)] + 0{ul)^ . 
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Here we have assumed in the second equahty that ouq becomes small in the vicinity of Tc and 
an expansion in cuq thus is appropriate. Similarly we can obtain the low frequency contri- 
bution to AGqq{t,T). Both correlation functions have the same functional dependence on 
Euclidean time for tT ^1/2. It is this region, where one can hope that the low frequency 
part of the spectral function dominates the r-dcpcndence of the correlation functions and 
thus will allow a determination of transport properties. 

Prom Eq. 31 we find that the low frequency part of the spectral function, which is propor- 
tional to bulk viscosity, contributes a term proportional to (^cuq to the value of the correlation 
function at the midpoint, tT = 1/2. As discussed in the introduction, universality argu- 
ments given for the scaling of the bulk viscosity in the vicinity of a critical point in the Ising 
universality class suggest, however, that the ratio of bulk viscosity and relaxation time, 
(/tr = (^cuq, vanishes as the inverse of the specific heat, i.e. it is expected to be proportional 
to the velocity of sound. In the vicinity of the correlation function Gqq{1/2T, T) thus is 
dominated by contributions from the large frequency, regular part of the spectral function. 
The sensitivity to bulk viscosity also is small away from the midpoint where contributions 
from the low frequency part of the spectral function are proportional to C^oJq. To separate 
this contribution from the large frequency contributions it would be necessary to become 
sensitive to the characteristic logarithmic dependence on Euclidean time that arises from 
the low frequency part of the spectral function. 

VI. CONCLUSIONS 

We have shown that the correlation function of the trace of the energy-momentum tensor, 
GeO) diverges at the critical point of the finite temperature SU(2) gauge theory in (3+1)- 
dimensions. The divergence is controlled by the critical exponent a of the 3-dimensional Ising 
model that also controls the divergence of the specific heat. Furthermore, we have shown 
that Gee becomes independent of Euclidean time at the critical point, which indicates that 
its spectral representation receives a contribution from a 5-function at zero frequency. This 
confirms the observations made in Ref. [10]. 

To get direct access to the bulk viscosity C and the characteristic frequency uoq that 
controls the relaxation time for density fluctuations will require much more detailed studies 
of the dependence of correlation functions on Euclidean time. In fact, we conclude that a 
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determination of bulk viscosity from correlation functions of the energy-momentum tensor 
is not possible without a simultaneous determination of the relevant relaxation time scale. 
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